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Abstract. - We study via renormalization group (RG) and large N methods the problem 
of continuum SU(N) quantum Heisenberg ferromagnets (QHF) coupled to gapless electrons. 
We establish the phase diagram of the dissipative problem and investigate the changes in the 
Curie temperature, magnetization, and magnetic correlation length due to dissipation and both 
thermal and quantum fluctuations. We show that the interplay between the topological term 
(Berry's phase) and dissipation leads to non-trivial effects for the finite temperature critical 
behavior. 



Introduction. - The subject of quantum ferromagnetism in metallic and semiconducting 
environments has recently attracted a lot of interest due to its relevance to a broad class 
of physical systems such as heavy fermion compounds, for instance URu2-£cRe :r Si2 51 anc ^ 
Thi_ :c U a ;Cu2Si2 [2], dilute magnetic semiconductors such as Gai-^Mn^As ferromagnetic 
dichalcogenides CeTe2 @], manganites Lai-^Sr^MnOs [5], and 2D electron systems in the 
quantum Hall regime |H| . These systems are characterized by the presence of local moments 
that couple directly to an itinerant electron liquid and can be described in terms of the 
Hamiltonian: 

U = - J J2 S i ■ S J + J K J2 Si ■ 4,0°™**,*' + Ke, (1) 

where J is a Heisenberg exchange between localized spins Si (total spin 5*), Jk is an ex- 
change coupling between localized spins and conduction electrons, c\ a (ci. a ) is the creation 
(anihilation) electron operator, a are the Pauli matrices, and TL e describes the conduction 
electrons. 

On the theoretical side, the critical properties of the insulating QHF have been studied 
via RG methods [7] and it has been shown that both the correlation length and the mag- 
netic susceptibility exhibit renormalized classical behavior at low temperatures and that no 
non-trivial quantum critical point (QCP) exists, in agreement with Monte Carlo data .$ . 
Other formulations of the insulating QHF based on O(N) and SU(N) models were used to 
explain magnetization and relaxation rate in incompressible quantum Hall systems |9lll)j . The 
problem of itinerant Hubbard ferromagnets, on the other hand, was studied by Hertz in 
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the context of classical Ginzburg-Landau- Wilson theory. This approach has been criticized 
recently because it does not take into account the physics of soft-modes |12| . 

In a metallic system such as the one described by Q dissipation occurs because electrons 
scatter from the localized moments and act as a heat bath for the spin dynamics. There 
are two main sources of dissipation: Landau damping as in the case of clean magnets, and 
electronic diffusion in the case of structurally disordered magnets |13| . We are going to use 
both momentum shell RG and large N analysis to generate the phase diagram of the dissipative 
problem. As we are going to show, a new QCP arises from dissipation. We find, moreover, that 
dissipation is a relevant perturbation under the RG and produces a runaway flow to strong 
dissipation. Our results from the RG analysis agree with the large N solution of the problem. 
Finally, we discuss the effects of dissipation on the finite temperature critical properties of 
the system Q. More specifically we discuss how the Curie temperature is affected by the 
different types of conduction electron dynamics. 

The effective field theory. - Our starting point is the partition function for the system 
described by Q (we use units such that h — ks = 1): 



is the long wavelength Euclidean action of the ferromagnetic O(N) non- linear sigma model |U] 
and f3 = 1/T. Here n is a N component vector that represents the local magnetization, 
p® = JS 2 a 2 ~ d is the spin-stiffness (a is the lattice spacing, and J is the nearest neighbor 
ferromagnetic exchange), and Mq = Sa~ d is the magnetization density in the ground state. 
The main difference between this action and the Hertz action, usually used for the study of 
itinerant quantum critical phenomena is the presence of the topological term (Berry's 
phase) described by the vector potential A of a Dirac monopole at the origin of spin space, 
a signature of the presence of localized moments. As we shall see, this term is fundamental 
for the understanding of the problem described by (JIJ. Furthermore, while Q describes a 
problem with short range interactions we can also introduce long range interactions that decay 
like 1 /r d+a by modifying the gradient term in momentum space from k 2 to k a with a < 2 |14| . 

The action Sj in J2J contains higher spatial derivative terms (with no time derivatives) 
corresponding to spin- wave scattering. It can be shown from a RG analysis that the coupling 
constant A associated with these terms has scaling dimension d + z — 4 and therefore is 
irrelevant for d < 4 — z, where z is the dynamical exponent UJ. Near the quantum critical 
point at T = we have z = 2 for an insulating magnet and z = 3(4) for a clean (dirty) magnet 
with local interactions |11|. This shows that A is irrelevant only if d < 2, 1,0, for z = 2,3,4 
respectively, otherwise producing corrections to the scaling of physical quantities. For the 
finite temperature T ^ critical phenomena, on the other hand, we can set z = and treat 
the temperature as a free parameter. In this case A is irrelevant in d < 4, and for this reason 
we ignore Si . We shall however keep z undetermined throughout our calculations in order to 
compare our results with other approaches. 



The third term in © contains the coupling to the dissipative environment and it has the 
form |13j : 




(2) 



where 






(4) 
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where uj n = 2irn/ (3, and 770 is the coupling to the electronic bath. The momentum dependence 
of the dissipative term depends on the electronic dynamics: Landau damping gives 5 = 1 (clean 
case) while diffusion implies 6 = 2 (dirty case) |13| . 

In what follows we rewrite the partition function (J2J) in terms of a basis of boson co- 
herent states (similar results can be obtained in the O(N) representation) ^U]. We write 
ni(r, x) = z q (t, x)cr^z / 3(r, x) where z a are complex bosonic fields subject to the condition 
that J2 a ZaZ a = N. In terms of these new fields So reads: 



where g = N/{2p° s ) is the bare coupling constant, A = iz a Vz a is a gauge field, and 
Co = 2Mogo/N is the bare topological constant (here the sum over a = 1,..,N is implic- 
itly assumed). Such formulation gives very good results at low temperatures even for N = 2 
but fails to reproduce the correct values for the critical exponents near T c |15|. However it is 
exact at all temperatures when N — > 00. Moreover, as we are going to show, our results for the 
effects of dissipation on the magnetism are already observed at the lowest temperatures, thus 
being robust, and should hold at higher temperatures as well, in particular at T c , irrespective 
of the type of formalism employed. 

While the non-interacting part of the action is quadratic in the z a fields, the dissipative 
term Sd turns out to be quartic. We can show, however, that in the large N limit the leading 
contribution to the boson self-energy is given by |lf>| : 



where (770) = i]o(zz). Thus dissipation in the scalar vector field n also induces dissipation in 
the bosons. 

Because of the constraint in J5J the problem described here is non-linear. Introducing 
a Dirac delta function for the constraint, the partition function can be written as Z sa 
J VtpT>zT>zT>Aevp |— S + J dr J d d x(f>(T, x)|, where (f>(r, x) is a Lagrange multiplier field, 
and 






(6) 



n.k. 





with 



Xo X ( k ) w n) = |k - A\ a - ic cj n + cj) + 7o-tJ-- (8) 



Here we have defined 70 = {^0)90/ N. Because of the quadratic form of S the fields z can be 
integrated out exactly and the partition function becomes Z = J £>02?_4exp {— NS e ff[(f>, A]} 
where 



9o Jo J 

When N — > 00 the partition function is dominated by the saddle point where (4>) = ro is 
constant and („4) = 0. The constant ro is determined by the saddle-point equation: 





n,k 



(10) 
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Notice that ro(T) = l/^ a (T) where £(T) is the magnetic correlation length. After summing 
over the frequencies we obtain: 

1 _ f A d d k [ n * dco k s 70 u n B (fiw) (n) 



9o J (27r) d 7_ o 7T {k s {k° + r -c uj)) 2 +76^2' 

where ns(z) = (e z — 1) is the Bose occupation number, A w 1/a is the high energy cut-off 
(for d = 3 we choose (Aa) 3 = 67r 2 for a cubic lattice so that the volume of the Brillouin zone 
is preserved), and — (A a /jo)k 5 is the bare frequency cutoff 

The RG calculation can be accomplished by splitting the momentum integral into two 
pieces: from < k < A/b and another from A/6 < k < A, with b = e e . After a change 
of variables, k = x 2 ' cr (A/b) and u = y(A/b) z , the parameters of the problem are rewritten 
as: r = r Q (A/b)-°, c = c (A/6) z — , 7 = l0 (A/b) z -^\ g = g Q S d (A/b) d+z - /a, and t = 
g S d (A/b) d - a /a/3, where Sj 1 = 2 d - 1 TT d / 2 T(d/2). In terms of the new variables, the saddle- 
point equation becomes: 

^ = 2 11{ X a ["" H ^^-'lynBjgy/t) 

where 

Z- 1 = 1 - - f d, . i4 _7^^ x /, (13) 

is a renormalization constant accounting for the contributions of the high energy modes. From 
Z we obtain the RG equations to leading order in 1/A 16 : 

di r = a r, 

di c = (a — z) c, 

dt 7 = [8 + a — z) 7, 

dig = -(d+ z - a) g + g 2 $(r, c, 7, 5, i), 

dit = -(d- ( 7)t + .gt$(r,c,7, 5 ,i), (14) 

with 

*(r,c™*) - - r ^ TTIT^TT ^ ' (15) 

i-f2 fc (1 + r - cy) 2 + ( 7 y) 2 

For <i > 4 — 2: we should add to the above set the equation di\ = (d + z — 4)A + p®/Mq which 
gives the flow of A [Hj ■ 

We hnd that there are three fixed points in the (t, g, 7) parameter space: (0, 0, 0) which is 
the trivial fixed point with r — (that is, with long range order); (f» = (d — a)/2, 0, 0) which is 
a unstable fixed point that marks the transition from ferromagnetism to paramagnetism; and 
a line of non-trivial fixed points at (0, g c (7), 7), with r = 0, that separates a fully polarized 
ferromagnet from a partially polarized ferromagnet |!Bv. Furthermore, we notice that the 
lower critical dimension (t» > 0) is given by di c = a. For short range interactions, a — 2, and 
order is only possible above two dimensions in accordance to the Hohenberg-Mermim- Wagner 
theorem. For long range interactions, a < 2 and long range order is possible for d < 2 because 
even long range distortions of the ferromagnetic order cost finite energy. 

We can now discuss the case of quantum critical phenomena when z ^ 0. For 7 — we 
recover the results for the insulating QHF by setting z = a so that c is unrenormalized under 
the RG while 7 is found to be a relevant perturbation. In fig. ^a) t c is the solution of 
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d(t = for 7 = and finite g, that is, it is the critical temperature of the insulating QHF 
and is depicted as the dashed line in fig.^b). For 7 ^ 0, on the other hand, we can choose 
z = 5+o so that 7 remains unrenormalized and we obtain z = 3(4) for the metallic ferromagnet 
with local interactions and Landau damping (diffusion). In this case, we see that c becomes 
irrelevant and flows to zero signaling the suppression of localized moments and making the 
renormalized action look very similar to Hertz action for the itinerant case. Although this 
result shows that c is irrelevant under the RG, one has to remember that the RG indicates 
the behavior of the system at very low temperatures, deep inside of the ferromagnetic phase if 
d > a. However, at finite energy scales (finite temperatures) c produces important changes in 
the behavior of T c as we explain below. Finally, we remark that a complete treatment of the 
RG for z ^ should also take into account the corrections from the flow of A to the scaling 
of all the other couplings [5]. 

In order to understand the critical behavior at finite temperatures we turn to the solution 
of saddle point equation . The critical line for ferromagnetic order is obtained by imposing 
that ro — 0, for j3 = (3 C - In fig. [5] we show the behavior of the Curie temperature as a function 
of the dissipation parameter rjo ■ We clearly see that the Curie temperature decreases faster for 
Landau damping than for diffusion. The physical interpretation is that ballistic electrons have 
an infinite mean free path and strongly scatter from the localized magnetic moments, thus 
introducing a large amount of fluctuations into the system. Diffusive electrons, on the other 
hand, have a finite mean free path and are consequently less effective in increasing fluctuations. 
We have also noticed an increase of about 1% on the value of the Curie temperature for 
small enough 770 1 with respect to the case without dissipation. The enhancement of the 
magnetism due to diffusion also seems to be important in the 2D MOSFET problem 17 . For 
Co — > we find that T c is rapidly suppressed by the dissipation as shown in fig. [21 showing the 
importance of the topological term for the finite temperature critical phenomena. Remarkably, 
as 7/0 — > 1, we observe an enhancement of the magnetism of about 20% with a T c higher in the 
diffusive (dirty) regime relative to the ballistic (clean) case. Real systems should then be found 
in between these two limits. We have calculated other physical properties of interest such as 
the magnetization, M , and the magnetic correlation length, £, as a function of temperature 
for the two limiting cases of Landau damping and diffusion for a fixed value of 770- The results 
are shown in figs.|3and0] As we see from fig.Olthc enhancement of the magnetism is already 
present at low enough temperatures where the formalism is more adequate. Consequently, 
such a behavior is expected to be observed within other formulations of the problem that 
captures the correct physics at high enough temperatures. 

In summary, we have studied the problem of metallic quantum ferromagnets at finite 
temperatures using large N and RG techniques and established the phase diagram for this 
problem. We found that from the point of view of the insulating ferromagnet, z — a, any small 
amount of dissipation produces a runaway RG flow to strong dissipation. We found a line 
of new QCP corresponding to a transition to a not fully polarized ferromagnet at g > g c , to 
be investigated elsewhere U^J. We also found that the topological and dissipative terms have 
a rather non-trivial competition that, as a consequence, produce a higher Curie temperature 
for diffusive than for ballistic electrons. Our theory provides a general and unified framework 
for the study of both insulating as well as metallic magnets (clean or structurally disordered) 
at vanishing and finite temperatures. Applications of our results to experimentally relevant 
systems will be given in future publications. 

* * * 

We acknowledge invaluable discussions with I. Affleck, K. Bedell, C. Chamon, P. Kopietz, 
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Fig. 1 - RG flow for d = 3, a — 2 and 8 = 2: (a) flow in the t x 7 plane; (b) flow in the t x g plane 
for g < g c . In (b) we compare the cases without (dashed line) and with (solid line) dissipation. 
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Fig. 3 - Reduced magnetization M/Mq as a function of temperature for S = 5/2 and r/o = 1.0. 
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Fig. 4 - Inverse correlation length as a function of temperature for S = 5/2 and r/o = 1.0. 



